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Abstract. We present evolutionary models of zero-metallicity very massive objects, with initial masses in the range 120 Mq 
- 1000 Mq, covering their quiescent evolution up to central carbon ignition. In the attempt of exploring the possible occurrence 
of mass loss by stellar winds, calculations are carried out with recently-developed formalisms for the mass-loss rates driven by 
radiation pressure (Kudritzki 2002) and stellar rotation (Maeder & Meynet 2000). The study completes the previous analysis by 
Marigo et al. (2001) on the constant-mass evolution of primordial stars. Our results indicate that radiation pressure (assuming 
a minimum metallicity Z = 10~* x Zq) is not an efficient driving force of mass loss, except for very massive stars with 
M ^ 750 Mq. On the other hand, stellar rotation might play a crucial role in triggering powerful stellar winds, once the 
nr-limit is approached. However, this critical condition of intense mass loss can be maintained just for short, as the loss of 
angular momentum due to mass ejection quickly leads to the spinning down of the star. As by-product to the present work, the 
wind chemical yields from massive zero-metallicity stars are presented. The helium and metal enrichments, and the resulting 
AY/AZ ratio are briefly discussed. 
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1. Introduction 

The first generation of stars ever born in the Universe 
(Population III) is assigned a role of paramount importance for 
several astrophysical issues. In particular, the formation of a 
primeval population of very massive objects (VMOs, with ini- 
tial masses in the range 10^ — 10^ Mq) was invoked in the '80s 
to explain several questions, such as the observed metallicities 
of Population II stars, the primordial helium abundance, the 
re-ionisation of cosmic matter after the Big-Bang, the missing 
mass in clusters of galaxies and galactic haloes, the G-dwarf 
problem (see CaiT et al. 1984 for an extensive review). In those 
years a few evolutionary calculations of zero-metallicity VMOs 
were carried out (e.g. Bond et al. 1984, El Eid et al. 1983, Ober 
et al. 1983, Klapp 1984ab), but afterwards the scientific pro- 
duction on the evolution of the first stars dropped off. 

The recently renewed interest in Pop-Ill stars, essentially 
driven by the nowadays flourishing of observations at very low 
metallicity and/or high redshift (refer to e.g. the proceedings 
of the ESO symposium edit by Weiss et al. (2000); see also 
Kudritzki et al. (2000), Bromm et al. (2001), Schaerer (2002), 
Panagia (2002) for extensive analyses on the expected observ- 
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able properties of primordial stellar populations) has again 
stimulated the calculation of stellar structures made evolve 
with initial metal-free chemical composition (e.g. Cassisi et al. 
2001; Marigo et al. 2001 and references therein). 

Marigo et al. (2001; thereinafter Paper I) presented an ex- 
tensive study of the evolutionary properties of zero-metallicity 
stars over a wide range of initial masses (0.7 Mq ^ M ^ 
IOOMq). In that work - to which the reader is referred for all 
the details - we calculated all stellar tracks at constant mass. 

Our study on Pop-Ill stars is now extended to very mas- 
sive objects (VMOs), with initial masses in the range 120 — 
1000 Mq, that is chosen not only for continuity with Paper I, 
but also in consideration of current theoretical indications that 
the primordial initial mass function (IMF) might have peaked 
in the (very) high-mass domain (with M ^ 100 M©; see e.g. 
Bromm et al. 1999, 2002; Nakamura & Umemura 2001; Abel 
et al. 2000). 

To this aim we calculate evolutionary models for zero- 
metallicity VMOs that cover the major core-burning phases, 
extending from hydrogen to carbon ignition. With aid of avail- 
able analytic formalisms, we address the question of the possi- 
ble occuiTence of mass loss via stellar winds, which is an im- 
portant, but still problematic, aspect of stellar evolution at zero 
metallicity. 
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In massive stars with "normal" chemical composition the 
principal driving force resides in the capability of metallic ions, 
present in the atmosphere, to absorb radiative momentum and 
transfer it to the outermost layers, that can be then accelerated 
beyond their escape velocity (Castor et al. 1975; Pauldrach et 
al. 1986). In a gas composed only of hydrogen and hehum, 
like the primordial one, the lack of metallic ions sets the first 
important difference and the question arises: Is the radiative ac- 
celeration due to the H and He lines strong enough to trigger 
significant mass loss ? Furthermore, may other possible pro- 
cesses - e.g. related to pulsation instability and stellar rotation 

- be efficient mechanisms to drive mass loss from Pop-Ill mas- 
sive stars ? 

Some of these questions have already been addressed 
by other investigations, others still deserve to be quantita- 
tively analysed and discussed. There are few studies in the 
literature that present evolutionary models of very massive 
zero-metallicity stars with mass loss during the pre-supernova 
phases. In the very massive domain (500 — 1000 Mq) Klapp 
(1983, 1984) carried out evolutionary calculations by adopting 
a simple empirical law for mass-loss (Barlow & Cohen 1977) 

- based on Galactic observations - that linearly scales with 
the stellar luminosity. On the basis of semi-analytical models 
of very massive objects (with masses 10^ — lO'^ Mq), Bond 
et al. (1984) argued that the same kind of dynamical insta- 
bility arising in the H-burning shell of Population I models 

- which should cause the ejection of the entire envelope - 
might also affect Population III VMOs, depending on the ac- 
tual abundance of the CNO catalysts in the H-shell. The evolu- 
tion of zero-metallicity stars (with masses 80 — 500 Mq) dur- 
ing the nuclear phases of H- and He-burning was calculated 
by El Eid et al. (1983) with the adoption of a semi-empirical 
relation to derive the mass loss rate (Chiosi 1981). Recently 
Baraffe et al. (2001, see also Heger et al. 2001) performed a lin- 
ear stability analysis on metal-free very massive models (with 
masses 120 — 500 Mq), and investigated the related possibility 
that mass loss may characterise the pulsation-unstable stages. 
Heger & Woosley (2002) calculated the evolution and nucle- 
osynthesis of quite massive Population III stars (with masses 
~ 100 — 300 Mq) including the final fate of the supernova 
event, but no mass-loss is assumed to occur during the hydro- 
static phases of major nuclear burnings. 

In the present work we attempt to further explore the issue 
of mass loss in primordial conditions, by carrying out new evo- 
lutionary calculations of zero-metallicity VMO with recently 
developed mass-loss formalisms, related to radiation pressure 
and stellar rotation. 

The paper is organised as follows. Section ||introduces new 
evolutionary calculations for zero-metallicity VMOs, having 
initial masses of 120, 250, 500, 750, and 1000 Mq. Mass loss 
is included with the adoption of recent analytic formalisms, 
namely: Kudritzki (2002) for radiation-driven mass loss, and 
Maeder & Meynet (2000) for the additional effect of stellar 
rotation. General characteristics of the models are analysed in 
relation to energetics, nuclear lifetimes, internal structure, sur- 
face and chemical properties. The efficiency of mass loss is 
discussed on the basis of the adopted formalisms. The corre- 
sponding predictions for the chemical yields, ejected via stellar 



winds, are given in Sect. ^ In particular, the expected helium 
and metal enrichment, produced by a hypothetical burst of Pop- 
III star formation, is derived and discussed as a function of the 
involved parameters. In Sect. ^ we express a few concluding re- 
marks. Finally, stellar isochrones for very young ages are pre- 
sented in Sect. 

2. Evolutionary calculations 

We have carried out evolutionary calculations of zero- 
metallicity very massive objects (VMO) with initial masses of 
120, 250, 500, 750, and 1000 Mq. All input physics are the 
same as those employed in Paper I - to which the reader is re- 
ferred for all details - except for the applied overshoot scheme. 
Overshooting is allowed from the borders of the convective 
cores with the parameter Ac = 0.5 (following Bressan et al. 
1981), whereas the classical Schwarzschild criterion is other- 
wise applied (i.e. convective shells and outer envelope). The 
evolution is followed starting from the zero-age main sequence 
(MS) up to the ignition of central carbon, i.e. covering the ma- 
jor stages of hydrogen and helium burning. 

In order to estimate the efficiency of mass-loss via stellar 
winds in primordial VMOs, we adopt recently-developed an- 
alytic formalisms (Kudritzki 2002; Maeder & Meyent 2000), 
whose basic ingredients are summarised in the next sections. 

2. 1. Radiation-driven mass-loss rates 

The role of radiation pressure - on the resonance lines of metal- 
lic ions like CIII, Nell, Oil - in driving mass loss from massive 
stars has been a firm result of theoretical astrophysics since the 
pioneer works by Lucy & Solomon (1970), and Castor et al. 
(1975; CAK). In short, the crucial term determining the effi- 
ciency of mass loss by line driven winds is the radiative line 
acceleration, which is expressed in units of the Thomson accel- 
eration via the so-called force multiplier. This latter, in turn, is 
suitably parametrised as Af [t(r)] = fci~"n*, where t{r) - the 
Thomson optical depth parameter - is proportional to the num- 
ber density of electrons, the Thomson cross-section, and the ra- 
tio between the local velocity of thermal motion to the spatial 
gradient of the stellar wind outflow velocity field; k represents 
the normalisation integral over the contribution of all photon 
absorbing spectral lines with different line strength; n is the 
ratio of the local number density of electrons to the local geo- 
metrical dilution factor of the radiation field; and finally a and 
5 are suitable exponents describing the optical depth and den- 
sity dependence of the radiative line force (see CAK, Abbott 
1982a, Pauldrach et al. 1986, Kudritzki, 1988, 1998, 2002 and 
references therein for more details). The quantities k, a, S are 
known as the force multiplier parameters. 

The expected dependence of radiative wind efficiency on 
metallicity (in the range from solar to 0.01 solar) has been 
investigated by e.g. Abbott (1982b), Kudritzki et al. (1987), 
Leitherer et al. (1992), Kudritzki & Puis (2000, and references 
therein), Vink et al. 2001. The mass-loss rates are predicted 
to decrease at lower metallicities as a power law following 
(Z/Z0)\ with 6 = 0.5-0.8. 
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Kudritzki (2002) in his investigation of winds, ionising 
fluxes and spectra of very massive stars at very low metaflic- 
ities has recently extended the analysis down to 0.0001 solar. 
These calculations introduce a new treatment of the radiative 
line force and are based on revised and updated atomic physics 
and line lists (see e.g. Puis et al. 2000 and references therein). 
Following a more realistic consideration of the line strength 
distribution function, the usual assumption of constant force- 
multiplier parameters (k, a, S) throughout the atmosphere is 
replaced by the introduction of an explicit dependence of these 
parameters on optical depth and density. This improved treat- 
ment of the line force is needed at extremely low metallicity, 
because the depth dependence of the radiative acceleration dif- 
fers significantly from the case of higher metallicity. As a con- 
sequence, a new elaborate algorithm to calculate the conditions 
at the critical point of the wind is also needed. 

The results of Kudritzki's stellar wind calculations are fit- 
ted by means of analytic formulas, that give the mass-loss rate 
as a function of effective temperature, luminosity and metal 
abundance for very massive stars with masses above 100 Af©: 



logM = gi([Z]-[Z,,i„])"-5 + Q„ 



inAfgyr-i (1) 



where we define 



Ql — Qo ^ Qmin(— [^]min) 



[Z] - log(Z/Z0) (2) 



with Z being the metallicity (abundance in mass fraction). The 
quantities [Z],nin, Qmin, and Qq are expressed as polynomials 
of the luminosity, in the form 

/ = ao + aiL + asL^ with L = log(L/L0) - 6.0. (3) 

The coefficients are given as a function of the effective tem- 
perature (see table 3 in Kudritki 2002), in the range from 40000 
to 60000 K. 

The above formalism is applied with the following assump- 
tions: 

- At any time the metallicity - that enters Eqs. (|lj)-(||) - is 
set Z — max{Zj-nin, Zctf)- This means that when eval- 
uating the radiative mass-loss rate we do not assume a 
truly zero metal content, but the maximum between the 
lowest value in the validity range of Kudritzki's formula 
(Znun = 10~* X Zq), and Z^s ^ 1 - X -Y. This lat- 
ter denotes the effective current metallicity, with X and Y 
being the hydrogen and helium abundances, respectively. 
In this way we may account, to a first approximation, for 
any possible surface enrichment in metals (mainly CNO el- 
ements, such that Zcff > Zj^in) exposed to the surface by 
stellar winds or convective mixing. 

- For Tos > 60000 K, we choose not to extrapolate the fit 
formula beyond its validity domain, but constrain T^s = 
min{Tca, 60000). It should be noticed that most of the 
main-sequence lifetime of VMOs is spent at Toff > 60000. 

2.2. Effect of rotation 

As shown by e.g. Friend & Abbott (1986) and Pauldrach et al. 
(1986), radiative mass-loss rates are enhanced by rotation by 



a factor that depends on the ratio, Wrot/i'crit, between the rota- 
tional velocity (usually taken at the equator), and the critical or 
break-up velocity. This latter is defined by the vanishing of the 
net radial acceleration contributed by the gravitational, radia- 
tive, and centrifugal forces. 

In our work we apply the analytic recipe presented in 
Maeder & Meynet (2000, and earlier works of the series refer- 
enced therein), where the results of a detailed quantitative anal- 
ysis of stellar rotation are conveniently synthesised. According 
to Maeder & Meynet's scheme, the mass-loss rate from a rotat- 
ing star can be expressed as a function of the purely radiative 
Mr ad, and the rotational correction factor: 



M(i;,ot) ^Fnx M,ad(frot = 0) , 
with the correction factor given by 



Fn 



(1-r)^ 



[i-(To + r)]- 



(4) 



(5) 



where F is the Eddington factor for electron screening opacity, 
is the angular velocity, pm is the internal mean density, and 
a is the force multiplier parameter corresponding to the slope 
of the line strength distribution. 

The effect of rotation is contained in the term 



2nGp„ 



4<i 
9 

^ '^crit 



where the break-up velocity is evaluated as 

2GM\3 



(6) 



(7) 



with clear meaning of the quantities. 

We refer to the classical il-limit when the centrifugal accel- 
eration balances the gravitational one, that is when t^rot = Vdit 
and the effects of radiation can be neglected. This is not the 
case of VMOs, as they are characterised by large Eddington 
factors. The so-called ilF-limit applies, instead, when the ef- 
fects of both radiation and rotation concur to balance the grav- 
itational acceleration. This corresponds mathematically to the 
divergency of the factor Fn, that is when the denominator of 
Eq. (^) vanishes. 

It is worth noticing that for Vi-ot = i^crit, the denominator 
of Eq. (g) becomes - [0.361 - F] and the 17F-limit could 
be met for values F > 0.639, that is for luminosities fainter 
than the Eddington limit. Viceversa, if F > 0.639, the f2F- 
limit could be already reached for rotational velocities lower 
than Wcrit, that is below the i7-limit. 

If f2F-limit is attained, then Eq. (^ becomes meaningless, 
and the theory does not indicate any straightforward simple 
treatment for such critical occurrence. For the sake of simplic- 
ity, we handle this circumstance by assuming a constant value 
for the mass-loss rate, Md-a, its value being a free parameter to 
be specified in our calculations. 

The evolution of the rotational velocity v^ot = ^{R) x R 
(where R is the surface radius) is derived under simple assump- 
tions, namely: 

- Spherical symmetry 



p. Marigo et al.: Evolution of zero-metallicity VMO with mass loss 



4 

- Rigid-body rotation, ri(r) = 

- Conservation of the current total angular momentum over 
any evolutionary time step 

- The equations of stellar structures are kept unchanged, i.e. 
without including any rotational term. 

In summary, we proceed as follows. Let us suppose that at a 
given time t the total angular momentum of the star is 

c{n) = nxj, (8) 

where 

J^^rp{r)r'dr (9) 
•J Jo 

is the moment of inertia of the stellar sphere. All quantities, like 
the stratification in radial coordinate r, density p and surface 
radius R, are evaluated at t. 

If during the evolutionary time step, St = t' — t, the star 
has lost a shell of mass SM, comprised between radii f and R, 
then the total angular momentum at time t' is given: 

c'{n') = n' X — r p'ir'y^dr' (lo) 

3 Jo 

o rR 

= C{n)^nx ^ / p{r)rMr 

where all primed quantities are now referred to time t'. 

In this way we can account for the decrease of the total an- 
gular momentum due to mass loss (see Heger & Langer 1998), 
as well as the effects of possible spinning-up (slowing-down) 
of rotation due to contraction (expansion) of the star 

The free rotation parameters to be specified before calcu- 
lating each model are: i) the initial rotational velocity v,-ot,o, 
and ii), the force multiplier parameter a, and the critical mass 
loss rate A/crit at the i7r-limit. In our calculations we adopt 
Wrot,o = 500 km s-i; a = 0.1, 0.4; and Afcrit = 10"^ Mq 
yr~^. 

We set a lower than the standard value of 2/3 valid for 
hydrogenic ions since, according to Kudritzki (2002), a is 
found to decrease at decreasing metallicities. We notice that the 
lower a, the larger the correction factor due to rotation Fn (see 
Eq. (^)). However, since our results for the case a = 0.1 do not 
substantially differ from those for a = 0.4, we will show and 
discuss only these latter. 

Finally, we choose to set Merit of the same order of mag- 
nitude as the very large mass-loss rates measured in Luminous 
Blue Variables (LBV), yellow and blue super-giants (see e.g. 
Leitherer 1997). 

2.3. General characteristics 

Let us now first consider the general features that characterise 
the evolution of zero-metallicity VMOs, and then move to dis- 
cuss the efficiency of mass loss. 

It is well known that in stars originally lacking in CNO ele- 
ments the (pre-)main sequence gravitational contraction of the 
central regions may proceed to such an extent that the density 
and temperature conditions allow the ignition of the triple— a 




10 100 1000 

Fig. 1. CNO and ^^C central abundances (in mass fraction) 
when 5% of hydrogen has been burnt in the convective core, 
as a function of the initial stellar mass 

reaction, with consequent activation of the CNO-cycle. In stars 
with M ^ 20 Mq, the first occurrence of the triple— a reac- 
tion usually takes place after core H-burning has already started 
via the p-p reactions, the time delay between the two events re- 
ducing at increasing stellar mass. Eventually, in more massive 
stars, with M ^ 20 Mq, the first synthesis of primary ^^C 
is practically simultaneous with the onset of core H-burning, 
which proceeds via the CNO cycle since the very initial stages. 

This is illustrated in Fig. |l], showing the primary CNO and 
^^C abundances in the convective core at the beginning of H- 
burning, as a function of the initial stellar mass. It is worth 
noticing that X(CNO) ~ X(^'*N), because as soon as the 
triple-a reaction produces some ^^C, most of it is immediately 
burnt into ^''N, given the high temperatures at which the CNO- 
cycle operates. For M > 20 Mq, the total X{CNO) has al- 
ready grown to ^ 10^^*^ — 10^^, which guarantees the primary 
role of the CNO-cycle as nuclear energy source for the star dur- 
ing the whole main sequence phase (see Paper I). The abrupt 
drop of X(CNO) at lower masses indicates that the full activa- 
tion of the CNO-cycle is delayed to later stages of H-burning, 
initially proceeding via the p-p reactions (see also figure 5 in 
Paper I). 

Various model predictions for relevant quantities are pre- 
sented in Table [l|. Let us first consider the nuclear lifetimes of 
core H- and He-burning (th and rne) as a function of the ini- 
tial stellar mass. In both cases the trend is initially decreasing 
at larger masses, until it flattens out towards a nearly constant 
value of roughly 2 x 10^ yr for th and 2 x 10^ yr for thc, 
towards the most massive models (M ^ 500Mq). It is worth 
briefly recalling the reason for this, with the aid of approximate 
scaling relations (see Kippenhahn & Weigert 1994), e.g. the 
definition of the nuclear H-burning lifetime, th oc M/ L, and 
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Table 1. Evolutionary properties of zero-metallicity VMOs for the adopted mass-loss prescriptions. From left to right the table 
entries read: the stellar initial mass; the H- and He-lifetimes; the mass of ejecta, carbon and helium cores at central C-ignition; 
the wind yields of helium, and C,N,0 elements. All masses are expressed in Mq, including the chemical yields. The nuclear 
lifetimes are given in yr. 



M 








Mc 




yC^He) 




yC^'N) 


yC'O) 


MODELS with v^ot = 


0; M = A/,ad 














120 


2.745E+06 


2.557E+05 


0.02 


50.94 


59.95 


4.090E-03 


6.805E-18 


5.536E-16 


1.056E-17 


250 


2.219E+06 


2.281E+05 


0.19 


117.29 


132.93 


8.655E-07 


2.956E-16 


1.966E-14 


3.595E-16 


500 


1.943E+06 


2.132E+05 


11.38 


243.00 


263.89 


1.725E-02 


5.251E-12 


4.870E-10 


8.745E-12 


750 


1.852E+06 


2.111E+05 


146.20 


343.84 


390.89 


1.324E+01 


4.229E-09 


3.983E-07 


4.842E-09 


1000 


1.861E+06 


2.153E+05 


581.29 


346.17 


418.67 


1.733E+02 


4.109E-01 


2.057E-05 


3.089E-01 


MODELS withurot.o 


= 500 km s" 


1 














120 


2.745E+06 


2.550E+05 


1.43 


50.92 


59.99 


1.587E-08 


3.039E-18 


2.478E-16 


4.728E-18 


250 


2.219E+06 


2.283E+05 


3.35 


117.50 


132.70 


3.853E-06 


1.085E-15 


8.520E-14 


1.587E-15 


500 


1.941E+06 


2.169E+05 


20.67 


245.22 


268.93 


5.848E-01 


1.767E-10 


1.507E-08 


2.877E-10 


750 


1.853E+06 


2.092E+05 


158.30 


343.63 


385.81 


1.612E+01 


5.002E-09 


4.682E-07 


5.742E-09 


1000 


1.862E+06 


2.155E+05 


581.41 


344.27 


414.40 


1.697E+02 


1.184E+00 


1.733E-05 


1.501E+00 



the mass-luminosity relation, L oc Af, so that th oc M^^^ . 
As the exponent 77 is typically 3-3.5 (say in the range from 
about 1 to 100 Mq), the nuclear lifetime is expected to decrease 
at increasing stellar mass. The almost invariance of th seen at 
very large masses is the result of the growing contribution of ra- 
diation to the total pressure, which makes the mass-luminosity 
relation approach a linear proportionality {rj ^ 1 and L cx M), 
so that Th ~ const. 

One common feature of the VMO models is the develop- 
ment of very large convective cores during both the H- and He- 
burning phase, with a typical fractional extension amounting 
to 70 - 90 % of the total mass. Moreover, at the exhaustion 
of hydrogen in the core central convection does not disappear, 
and the onset of core He-burning immediately follows because 
of the already high temperatures attained (logT ~ 8.2 — 8.3). 
This implies that the intermediate phase of H-shell burning is 
not expected for these stars. 

Two examples are illustrated in Figs. ^ and ^ that show 
the evolution of both convective and burning regions in two 
models with initial masses of 250 Af© and 1000 Af0, for two 
choices of the mass-loss prescription. From these plots we can 
already see that new products synthesised by nuclear burnings 
may be exposed to the stellar surface because of i) the inward 
penetration of the convective envelope across a chemical profile 
left by the gradual recession of the convective core, and/or ii) 
the progressive peeling of the stars due to mass loss. This point 
will be discussed in more detail in the following sections. 

Figure ^ displays the evolutionary tracks in the H-R dia- 
gram for the two sets of VMOs, calculated adopting the pre- 
scriptions for purely radiative and rotational-radiative mass 
loss, respectively. In all cases central H-burning is mostly spent 
at high effective temperatures, regardless the efficiency of mass 
loss. We can notice that the locus of points describing the on- 
set of core He-burning (practically coincident with the end of 
core H-burning) bends towards lower Toff at increasing stel- 
lar mass. On the other hand, the He-burning phase takes place 



in different regions of the H-R diagram, depending on stellar 
mass and mass loss. Soon after the onset of He-burning, models 
with I2OA/0 < M < 750 are able to reach their Hayashi line, 
where they remain up to central carbon ignition. The evolution 
of models with A/ = 1000 Mq takes place entirely at high ef- 
fective temperatures, essentially due significant mass loss. Next 
section will examine this aspect in mode detail. 

Finally, Fig. |^ shows the expected behaviour of the mass of 
the He- and CO- core at the stage of central carbon ignition, as a 
function of the initial mass in the domains of massive and very 
massive stars (8 Mq < M < 1000 Mq). The relation is almost 
linear, without any important effect due to mass loss. Our pre- 
dictions indicate that the interval A/ho = 64 to 133 Mq, which 
would produce pair-instability supernovae according to Heger 
& Woosley (2002)'s recent calculations, correspond to stellar 
progenitors with initial masses M = 127 to 252 Mq (see also 
El Eid et al. 1983, Ober et al. 1983, Bond et al. 1984, Carr et 
al. 1984). 

2.4. Mass-loss efficiency 

We will discuss now the efficiency of mass loss in relation to 
the two driving mechanisms here considered, namely: radiation 
(Sect. 2.1) and rotation (Sect. 22). Stellar ejecta are presented 
in Tabled 

Models with M = Afiad- In general, the purely-radiative 
mass-loss rates remain quite low during the whole evolution of 
VMO models with M < 500 A^©. For instance, the 120 Mq 
model returns only 0.02 Mq to the interstellar medium. In this 
mass range, the general properties are practically the same ones 
as those predicted for constant-mass evolution. As reasonable 
extrapolation we could expect that at lower masses (< 120 Mq, 
hence lower luminosities), the inefficiency of mass loss should 
even more marked, which supports the scenario of constant- 
mass evolution of Pop-Ill stars discussed in Paper I. Another 
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Fig. 2. Evolutionary properties of 
the 250 Mq model with the pre- 
scription for the purely radiative 
mass loss. Top panel: Stellar lu- 
minosity and effective temperature 
as a function of time during the 
major nuclear burnings up to cen- 
tral carbon ignition. Bottom panel: 
Convective and burning regions (in 
mass coordinate from the centre to 
the surface) across the stellar struc- 
ture. The upper solid line represents 
the mass coordinate of the stellar 
surface 
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Fig. 3. The same as in Fig.g, but 
referring to the evolution of the 
1000 Mq model calculated with the 
prescription for radiative-rotational 
mass-loss. The assumed initial ro- 
tational velocity is v^ot — 500 
km s^^. An almost identical plot 
is obtained for the model with the 
same initial mass, but assuming the 
purely radiative mass-loss rates 



common characteristic of the models with M = 120 — 500 Mq 
is the morphology of the tracks in the H-R diagram, and hence 
the location of major nuclear burnings. Actually, if the H- 
burning phase is spent in regions of high effective temperature 
(typically at log Toff ^ 4.8 — 5.0), the subsequent He-burning 
phase takes place far away from there, as the star evolves along 
its Hayashi line (typically at logTcff ~ 3.7 — 3.8). Indeed, the 
tracks of models with M < 500 Mq are characterised by quite 
a large excursion in Toff, that is performed at almost constant 
luminosity. 

Similar considerations on the morphology of the tracks in 
the H-R diagram apply to the 750 Mq model as well, though 



in this case the efficiency of mass loss is somewhat more sig- 
nificant, the total ejecta amounting to about 20 % of the initial 
mass (see Table [l|). 

A different behaviour characterises the most massive model 
here considered, with M — 1000 A/©. In this case, radiative 
mass loss is predicted to be strong since the very beginning 
of the MS phase, due to the combination of large luminosities 
and effective temperatures. The entire evolution is confined in 
regions of the H-R diagram of high effective temperatures, and 
the corresponding stellar ejecta up to the central carbon ignition 
are considerable, reaching to more than half the initial mass. 
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Fig. 4. Zero-metallicity stellar tracks of VMO with initial 
masses of 120, 250, 500, 750, 1000 Mq. The evolution covers 
the H- and He-burning phases up to ignition of central carbon. 
Top panel: Purely radiative mass-loss prescription; Bottom 
panel: Rotational-radiative mass-loss prescription 



The net result is that radiation-driven stellar winds do not 
produce appreciable mass loss from massive zero-metallicity 
stars, except for the most massive models, with M > 750 Mq. 
Our finding supports previous indications in this sense by 
Kudritzki (2000) and Bromm et al. (2001). 

It is worth making a final remark by looking at Fig. |], which 
displays the degree of chemical self-pollution at the surface, in 
terms of CNO abundance, produced by the combined effect of 
mass loss and convective mixing. As we see in most models the 
maximum enrichment does not exceed 10^^ — 10^^, that is the 
critical CNO abundance required in stars with original metal- 
free composition to fully activate the CNO-cycle as dominant 
energy source (see Sect. 13 and Fig. [l|). This occurs when the 
mass-loss front is only able to penetrate into the chemical pro- 
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Fig. 5. Surface abundance of CNO elements (in mass fraction) 
as a function of the evolutionary time, from the ZAMS up to 
carbon ignition in the core. Curves correspond to different ini- 
tial masses (as indicated) of VMOs with initial zero metallic- 
ity. The radiative mass-loss prescription is applied. Note that in 
most cases X(CNO)surf remains lower than Zmin assumed in 



the radiative mass-loss formula (see Sect. 2. 1 ) 



file generated by the recession of the convective core during 
the MS phase, as it is the case of the 500 and 750 Mq models. 
Models with lower masses, 120 and 250 Mq, show compara- 
ble self-enrichment, that is however reached after the end of 
the He-burning phase, not as a consequence of mass loss but 
because of deep convective dredge-up while on the Hayashi 
line. It follows that such degree of surface chemical enrich- 
ment should not produce enough resonant metal Unes in the 
atmosphere to activate efficient radiation-driven mass-loss. 

On the contrary, the most massive 1000 A/0 model exhibits 
a remarkable increase of the CNO abundance during the He- 
burning phase, due to the progressive reduction of the stellar 
mass. As shown in Fig. |] the H-burning shell is extinguished 
and eventually, close to the exhaustion of central helium, the 
mass-loss front penetrates into regions previously reached by 
convective core. Then, significant amounts of carbon and oxy- 
gen are exposed at the surface (X(CNO)« 0.6). Whenever this 
circumstance may trigger substantial radiative mass-loss is an 
interesting point to be considered, and it would deserve a dedi- 
cated study in the framework of the radiation-driven-winds the- 
ory. Finally, it is worth remarking that an additional chemical 
pollution of the stellar surface may be caused by rotationally- 
induced mixing (e.g. Heger & Langer 2000; Meynet & Maeder 
2002), which is not taken into account in the present analysis. 

Models with M = M{viot)- The amplification effect on mass 
loss due to stellar rotation turns out rather limited. As we can 
see, the evolutionary tracks in bottom panel of Fig. ^are almost 
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Fig. 6. Evolution of rotational parameters referring to the 120 
Mq model, calculated with T^ot.o = 500 km s~^, and A/crit = 
10~^ Mq yr^^. Top-left panel: Relevant terms determining the 
rotation correction factor given by Eq. The il- and VLT- 
limits are indicated by horizontal lines delimiting the maxi- 
mum values for Tn and + F, respectively. Top-right panel: 
Rotation correction factor. Bottom-left panel: Evolution of sur- 
face (Wrot) and critical (wcrit) rotation velocities. Bottom-right 
panel: mass-loss rate 

indistinguishable from those in the top panel. Actually, the to- 
tal ejecta are larger than those for models with radiative rates, 
but not such to drastically alter the evolutionary properties just 
presented. 

To better understand the results we refer to Figs. ^ and ^ 
showing the evolution of rotational parameters and mass-loss 
rates. In both cases the correction factor Fq (Eqs. (Q) and (||)) 
keeps close to one during nearly the whole MS phase, so that 
the mass-loss rate remains quite low (s^ 10^* Mq yr~^). 

Eventually, during the stages just preceding and following 
the H-consumption in the core, the 120 A/© model reaches both 
the ri- and ilF-limits at the same time. This result is supported 
by recent analyses of rotating stars at low metallicities. In fact, 
as discussed by Meynet & Maeder (2002) the approach of the 
i7F-limit is just favoured in massive and luminous stars at very 
low- or even zero metallicity. This circumstance should be the 
consequence of two concurring factors, namely: i) the nearly 
total conservation of the initial angular momentum during the 
almost constant-mass evolutionary stages, that would lead to 
the natural settling of the break-up rotation, hence the 51-limit; 
and ii) the large values of the Eddington factor in massive and 
luminous stars. In practice, both factors favour the divergency 
of the correction factor given by Eq. (Q). 

With respect to the former point i), we recall that various 
analyses (e.g. Heger & Langer 1998, 2000; Meyner & Maeder 
2000, Maeder & Meynet 2001) have shown that a massive 



Fig. 7. The same as in Fig. |, but for the 500 Mq model 



star evolving at (nearly) constant mass could approach or even 
reach the break-up rotation, towards the end of the MS phase 
(or in a blue-loop off the Hayashi line). In fact during these 
stages, the progressive stellar contraction would result in in- 
creasing surface rotational velocities (the so-called "spin-up 
phases"). This is exactly what we find in our 120 Mq model 
displayed in Fig. |[ At this stage the stellar mass is practically 
the same (^ 119.98 Mq) as the initial one. 

However this critical regime, that we describe by assuming 
an extremely high mass-loss rate Mait, is maintained only for 
very short. In fact, the reduction of total angular momentum 
carried away by the outermost ejected layers, and the increase 
of the radius due to the subsequent expansion (following the 
beginning of core He-burning) make rotation to slow down (see 
bottom-left panel of Fig. ||), with consequent departure from 
the Q. and riF-limits. 

Then, the star evolves rapidly towards lower effective tem- 
peratures onto its Hayashi line, which causes a further drastic 
decrease of the rotational velocity. The total ejecta with the in- 
clusion of the rotational effect is almost negligible, ^ l.SM© 
(see Table |l]). 

A similar behaviour is found for the 500 Mq model. 
However, in this case the larger luminosity concurs to lead the 
star closer to the ilF-limit already during the MS phase. Then, 
the radiative mass-loss rate is enhanced by a few times (see top- 
right panel of Fig. Again, the ilF-limit (but not the r^-limit) 
is briefly touched during the spin-up phase towards the end of 
the MS phase, but soon after abandoned because of the loss of 
angular momentum and the expansion of the star. Compared to 
the model with the same initial mass but with radiative mass 
loss, the net effect of rotation is to increase the total ejecta by 
roughly a factor of two, which still remains very small (see 
Table |). 
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At larger masses {M = 750, 1000 A/q) the effect of ro- 
tation turns out even weaker In fact, for these models the 
radiation-driven mass loss becomes efficient, which leads to the 
removal of a significant part of the original angular momentum 
already during the MS, hence preventing the approaching of 
the critical Vt and Q, — F-limits. 

2.5. Cautionary remarks 

At this point it is worth making a few cautionary remarks about 
the assumptions and simplifications adopted in our work to de- 
scribe the effect of mass loss in primordial massive stars. 



As far as the radiative mass-loss is concerned (Sect. 2.1 ), it 
should be pointed out that the position Z = max {Z nun, ^cff) 
and the treatment for T^b > 60000 K may lead to over- 
estimate the actual mass-loss rates, since for instance the min- 
imum metallicity, Zmm ~ 10~^, is typically two orders of 
magnitude larger than ZcS, (except for the 1000 Mq model). 
Moreover, the application of the formula for ilfrad to mod- 
els with M > 500 A/q requires to extrapolate the predictions 
above the upper-limit in luminosity (i.e. log(L/LQ) = 7.03) 
of the validity domain. 

The simplifying assumption at this stage, without calcula- 
tions for Teft > 60000 K and log(L/L0) > 7.03 available, is 
that the line force parametrisation at higher temperatures and 
luminosities is still reasonably approximated in this way. New 
calculations extending the temperature and luminosity domains 
are certainly needed for the future. 



As far as the rotational mass-loss is concerned (Sect. 2.2), 



we recall that the results are obtained in the ambit of a very sim- 
plified description of stellar rotation (mainly due to the assump- 
tion of solid-body rotation). It is clear that a strictly correct pro- 
cedure should consider important processes like meridional cir- 
culation, shears, and horizontal turbulence, which would also 
determine a feed-back on stellar structure (see e.g. Kippenhahn 
& Thomas 1970; Meynet & Maeder 1997). Among the main 
related effects, it turns out that the outward transportation of 
angular momentum from the central regions by meridional cir- 
culation would produce a differential rotation. These aspects 
are not included in the present work. 

Rather, our assumption of solid-body rotation would ideally 
correspond to the case of maximum efficiency of convection 
and shears in the transport of angular momentum from the inner 
to the outer regions of the star In this way even a relatively 
small loss of material from the surface leads to a substantial 
loss of the total angular momentum, and consequently to the 
slow-down of rotation. 

On the other hand, in the context of a detailed treatment 
of rotation - coupled to stellar structure and including also 
the role of meridional circulation - we may expect that the 
outward transport of angular momentum is less efficient. In 
fact, detailed computations accounting for differential rotation 
(e.g. Meynet & Maeder 2000, 2002) indicate that i) the angu- 
lar velocity increases progressively with time towards the cen- 
tre, and ii) the break-up rotation could be reached earlier and 
maintained for longer periods. Moreover, the possible effect 
of anisotropic mass loss (Maeder 2002) could further reduce 



the loss of angular momentum and lead to earlier break-up. 
Finally, two additional effects are expected to gain increasing 
importance at lower metallicity, namely: i) faster core rotation 
and more efficient rotational mixing (Meynet & Maeder 2002), 
and ii) possibly larger initial rotational velocities (Maeder et al. 
1999). 

These considerations lead to the possibility that the com- 
bined effects of rotation and mass loss may be somewhat under- 
estimated in our calculations. 

It follows that the inclusion of all these aspects might pro- 
duce different results compared to those presented in this pa- 
per. However, since up to now there is no dedicated analysis 
addressing such issues specifically for Pop-Ill stars, we believe 
that our simplified exploration of the effect of rotation on pri- 
mordial stars should keep its own validity. 

3. The stellar yields 

If ever formed, massive primordial stars marked the first step 
in the chemical evolution of the Universe (e.g. Oh et al. 2001; 
Scannapieco et al. 2002). The efficiency of their contribution to 
the chemical enrichment of the pristine gas depends essentially 
on two factors, namely: the individual stellar yields, and the 
details of the IMF. 

As for the former ingredients, predictions of elemental 
yields are available in Woosley & Weaver (1995) for super- 
novae II explosions produced by zero-metallicity stars with ini- 
tial masses in the range 11 — 40 A/q, and in Heger & Woosley 
(2002) for pair-instability supernovae, involving helium core 
masses in the range 64 — 133 Mq (corresponding to initial stel- 
lar masses of 140 — 260 Mq under the hypothesis of constant- 
mass pre-supernova evolution). According to the latter work, 
at both lower and larger masses, no explosive nucleosynthe- 
sis should be injected into the primordial interstellar medium 
(ISM) but, rather, trapped into black holes. For past works on 
pair-instability supernovae the reader could refer to El Eid et 
al. (1983), Ober et al. (1983), Bond et al. (1984), Carr et al. 
(1984). 

In their evolutionary calculations at initial zero metallic- 
ity Woosley & Weaver (1995) and Heger & Woosley (2002) 
did not consider any mass loss during the hydrostatic phases 
(except for that driven by pair-instability pulsation). Past cal- 
culations of the wind contributions from VMOs can be found 
in the works by El Eid et al. (1982), Klapp (1983,1984), and 
more recently Portinari et al. (1998). These latter estimates 
are based on the results of semi-analytic stellar models devel- 
oped by Bond et al. (1984). Our predictions of the wind ejecta 
from Pop-Ill VMOs are given in Table |l] that presents, for each 
model with initial mass Af , the net wind yields (in A/q) of ^He, 
i^C, ^"^N, and ^''O, the total ejecta AA/cj, and the masses of 
the He-core (A^hc) and C-core (Mc) at the ignition of central 
carbon. For the sake of completeness, the results are grouped 
on the basis of the adopted mass-loss prescriptions, namely: 
purely -radiative mass loss, and radiative-rotational mass-loss. 
However, the differences between the two sets of yields are 
quite small (see Sect. 2.4 ). 

It turns out that the wind yields ejected by VMOs are prac- 
tically negligible for all masses, except for the two most mas- 
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Fig. 8. Masses of the He- and CO- cores at the onset of car- 
bon ignition in the core, as a function of the stellar initial 
mass. They correspond to evolutionary calculations at con- 
stant mass for M < 100 M© (Marigo et al. 2001), and with 
radiation-driven mass loss for VMOs for 120 Mq < M < 
1000, AIq. Relevant ranges of A/hc are marked and labelled 
- following Woosley & Weaver (1995); Heger & Woosley 
(2002) - as a function of the final fate of the stellar pro- 
genitor, namely: i) black hole (BH) collapse without explo- 
sion (A^Ho > I33A/0); ii) pair-instability supernovae (PISN, 
64 Mq < A/ho < 133A/o); iii) black hole (BH) collapse after 
the final explosion (9 Mq < A/hc < 64Afo) iv) SNll explo- 
sion with stellar remnant (1.3 Mq < A^/hc < SAf©) 



sive models here considered, with M = 750, 1000 A^q. For 
instance, the 1000 Mq model with Wiot,o = 500 km s~^ ejects 
about 17 % of its initial mass as newly synthesized helium, and 
1 — 1.5 Af© in form of carbon and oxygen. In this particular 
case (see Fig. ||), during the MS phase the mass-loss front eats 
up the chemical profile left by the receding convective core and 
the wind ejecta, returned to the ISM, are mainly enriched in 
^He and ^^N. Then, during the subsequent He-burning phase, 
the deep peeling keeps going on to such an extent that, af- 
ter the ejection of the entire hydrogenic envelope (when the 
abundance of central helium is ^ 0.03), the mass-loss front 
extends beyond the extinguished H-burning shell, eventually 
exposing to the surface the nuclear products of partial core He- 
burning, mainly ^^C and ^^O. However, the net yields of these 
elements are not extremely large, given the short evolutionary 
time-scales involved from that stage up to the central carbon 
ignition (i.e. the end of calculations). 



3. 1. Chemical enrichment and primordial AY/AZ 

The extent of helium enrichment, AY, supplied by Pop- 
111 stars in the very early stages of cosmic chemical evolu- 
tion, is a potentially important issue for its possible cosmo- 
logic implications. For instance, the age of Globular Clusters 
(GC), as inferred from observations, depends sensitively on 
the assumed initial helium abundance Yp, which is commonly 
thought to sample the primordial Big Bang nucleosynthesis 
(BBN). Typically, an increment of the fractional mass abun- 
dance AFp = 0.02 would produce an age decrement of roughly 
15%, that is of the order of 2 Gyr within the range 10 — 15 Gyr 
(see e.g. Shi 1995). Moreover, the relative enrichment of he- 
lium with respect to that of metals - the so-called AY/AZ - 
produced by Pop-Ill stars, could largely differ from that derived 

- both theoretically and observationally - at non-zero metallic- 
ities, so that the usually adopted linear extrapolation down to 
Z = to infer the primordial helium abundance may lead to 
misleading results. 

Another cosmological issue is related to the nature itself of 
the dark matter present in galactic haloes and clusters. The stan- 
dard hot BBN predicts an upper limit to the baryonic density 
parameter, 51bbn < 0.06h^Q (where /150 is the Hubble con- 
stant in units of 50 km s^^ Mpc^^; see e.g. Kernan & Krauss 
1994), that is already lower than the present density param- 
eter inferred from dynamical estimates of virialized systems, 
^^dyn 0.1 — 0.2 (e.g. Dekel 1994). This discrepancy could be 
solved - without losing, at the same time, the consistency be- 
tween BBN predictions and the observed elemental abundances 
(like those of deuterium and lithium) at the lowest metallicities 

- by invoking that the bulk of dark matter is in a non-baryonic, 
still unknown, form. This is the commonly accepted picture. 

An alternative scenario might arise from the assumption 
that, instead, dark matter was made of baryons, e.g. in form of 
the remnants of pre-galactic stars. Were this the case, the im- 
mediate dramatic consequence would be that we should aban- 
don the standard BBN, and ascribe the first synthesis of light 
elements, like deuterium and helium, as well as the photon mi- 
crowave background to another source. Bond et al. (1983) and 
Carr et al. (1984) pointed at Pop-Ill very massive stars as po- 
tential candidates, and investigated the extreme possibility that 
these stars synthesised a substantial fraction, or even the total- 
ity, of primordial helium. In that case, a related constraint is that 
the accompanying metal enrichment does not exceed the upper 
limits set e.g. by the metallicities measured in Population 1 and 
11 stars. 

In order to derive an updated picture of the aforementioned 
aspects on the basis of the present results, we have calculated 
the helium and metal enrichments produced by a population 
of (very) massive zero-metallicity stars, with 10 Mq < M < 
8 X 10^ A^0. No other sources are considered, such as the wind 
yields from low- and intermediate-mass stars in the Asymptotic 
Giant Branch (AGB) phase (with 0.7 A/© <, M ^ 8Mq). In 
other words, we assume that the primordial initial mass func- 
tion favoured the formation of massive objects, a hypothesis 
that is supported by current models of primeval cloud fragmen- 
tation (e.g. Bromm et al. 1999, 2002; Nakamura & Umemura 
2001; Abel et al. 2000). The choice of 8 x lO"* Af© as maxi- 
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mum mass value is motivated by the indications that above this 
limit stars are not expected to experience static nuclear burn- 
ings because of post-Newtonian instabilities, and the chemical 
enrichment from these super-massive objects should be small 
(Fricke 1973, Ober 1979; see also Ober et al. 1983). 

In summary, we account for the chemical contributions 
from three groups of stars, namely: i) massive objects (MOs) 
with masses of 12-40 A/q, that experience no steady mass 
loss and undergo SN II explosions; and ii) VMOs with ini- 
tial masses of 120-260 Mq, that may experience mass loss 
via stellar winds and pair-instability SN explosions; and iii) 
VMOs with initial masses of 260 - 8 x 10"* M© that may suffer 
mass loss during the hydrostatic phases, before collapsing into 
black holes. As for the MOs, we adopt the supernova yields 
at Z calculated by Woosley & Weaver (1995). For the 
VMOs with initial masses of 120-1000 A/©, we combine our 
estimated wind contributions with the supernova yields (if non 
zero) by Heger & Woosley (2002). The matching between the 
wind and SN yields is performed as a function of the helium 
core mass as tabulated in Heger & Woosley (2002). Finally, for 
stars more massive than 1000 A/q, we adopt simple analytical 
prescriptions to estimate the helium and metal yields ejected by 
steady mass loss over the hydrostatic hydrogen-burning phase. 
Following Bond et al. (1984) we assume the maximum possi- 
ble helium yield, which is expressed as a function of the initial 
electron molecular weight (their equation 27). The correspond- 
ing metal enrichment is obtained by multiplying an average 
CNO abundance (newly synthesised during the core hydrogen- 
burning, typically (X(CNO)) « 10-9-10-8, see Sect. |]l|) by 
the amount of ejected mass (derived by equating their equations 
24 and 27). With our choices for the primordial hydrogen and 
helium abundances (Xq — 0.77, and Iq = 0.23 respectively), 
we obtain the fractional yields - normalised to the initial stellar 
mass - for objects with M > 1000 Mq: 
z/('*He)/Af = 0.17 
y(CNO)/Ar = 0.435 x(X(CNO)) 
In our calculations we assume (X(CNO)) lO^^. 

In a similar manner as in Ober et al. (1983), a very simple 
chemical model is assumed: stars formed in a single burst from 
a primeval cloud, that is assimilated to a closed box system, 
chemical enrichment is considered instantaneous with the burst 
and homogenised throughout the gas. 

The basic free parameters that define our model are: the 
fraction of the primordial cloud, with original mass Mq, glob- 
ally consumed in the burst of star formation, fs = M^/Mq, and 
the details of primordial IMF (in number of stars per unit mass 
interval): 



(/)(Af) = AM-"" 



for A/low < M < A/up , 



(11) 



where x is the slope of the assumed power law (x — 2.35 de- 
scribes the classical Salpeter relation); Afiow, M^p correspond 
to the minimum and maximum values of the initial stellar mass, 
respectively. 

The quantity A in Eq. ([ll|) must satisfy the condition 



Ms = / M(t){M)dM = A 



(12) 



Denoting with A^(A/), A^^(7\f ) the masses of the ejecta con- 
tributed by a star with initial mass M via stellar winds and 
supernova explosions, respectively, then the final mass of the 
cloud, left in form of gas after the burst of star formation, is 

/•A/up 

Mi = Mo - Ms + A [A^{M) + A^^{M)]M-''dM 

^ Mo-Ms + Ax n (13) 

Analogously, the final mass of the cloud in form of any given 
element j can derived from 

M, = {Mo - Ms)X,,o + A / {vYiM) + yf{M) (14) 

+ [A^{M) + A^^{M)]Xj^o}M^''dM 
= {Mo - Ms)Xj^o + Axnj 

where Xj o denotes the initial abundance of the element j, 
yY {M) and y|^(Af) are the corresponding yields (in mass 
units) due to stellar winds and SN explosion, respectively. In 
Eqs. ( [l3| ) and ( [l^ the quantities R and Rj are a measure of 
the integrated efficiency of the primordial stellar population 
in terms of restitution to the ISM of gas, and elemental abun- 
dances, respectively. 

Finally, after simple algebraic passages, we estimate the 
fractional abundance (by mass) of the element j in the gas re- 
sulting from the chemical enrichment by Pop-Ill stars: 



AXi = X, 



X 



_M^_ _ /s(7^J - nx. 



(15) 



We apply the above formula to derive the global enrichment of 
helium, AY, and metals AZ produced by a Pop-Ill burst of star 
formation. We set Yo — 0.23 and Zq = 0. The results can be 
analysed as a function of the efficiency of star formation (ex- 
pressed by the parameter fs), and different choices of the IMF. 
Figures ^ and |l0| show few illustrative examples, that bracket 
the widest variability range of the results for different combi- 
nations of the model parameters. In practice we consider five 
cases, namely: i)-ii) classical Salpeter IMF (with x ~ 2.35) 
over the mass intervals 10 - 1000 Af© and 10 - 8 x lO'* Af©; 
iii)-iv) Salpeter-like but flatter IMF (with x = 1.5) with the in- 
clusion of VMOs only, over the mass intervals 120 — 1000 Af© 
and 120 — 8 x 10'' Af©; and v) the same as the previous cases 
but for a very narrow mass range at about 1000 Af©. 

As for the helium enrichment, illustrated in Fig. ^ we can 
make the following remarks. In all cases, AY increases with 
the efficiency of star formation such that an amount of helium, 
comparable with the primordial value predicted by standard 
BBN, could be produced by Pop-Ill massive stars if /s ~ 1. 
This finding is in agreement with the results by Bond et al. 
(1983). However, leaving aside this extreme and unlikely con- 
dition, we notice in Fig. ^ that, at any fs, the largest contri- 
butions come under the assumption of an extremely top-heavy 
IMF either narrowly peaked at ~ 1000 Af© (short-dashed line), 
or including the contribution of super-massive objects up to 
8 X 10'* Af© (long-dashed fine). 

In these particular cases a helium enrichment already of the 
order of AY 0.01, hence able to affect the age determination 
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Fig. 9. Helium enrichment produced by a primordial simple 
stellar population as a function of the star formation efficiency, 
according to different prescriptions, as indicated 

of GCs, may be attained for /s w 0.05 — 0.1. Interestingly, these 
are typical values of the star formation efficiency as indicated 
by hydrodynamical models of structure and galaxy formation 
(e.g. Abel et al. 1999; Kawata 2001; Chiosi & Carraro 2002), 




Fig. 10. Helium enrichment as a function of metal production. 
The illustrated cases refer to the same prescriptions (and cor- 
responding line styles) as described in Fig. ^. Selected values 
of the star formation efficiency parameter, /s, are marked along 
the curves 



and chemical evolution of the intra-cluster medium (Moretti et 
al. 2002). Although the efficiency of star formation somewhat 
depends on the assumed IMF regulating the energy input from 
supernova explosions, stellar winds, UV flux etc, in reality it is 
the result of a delicate balance among several concurring phys- 
ical processes such as the energy input above, cooUng of gas 
by radiative mechanisms, and dynamical conditions establish- 
ing when a gas particle of the hydrodynamical models is able 
to form stars, so that the role of the IMF is of minor relevance 
(see Chiosi & Carraro (2002) for a detailed discussion of these 
topics). An efficiency of star formation /g « 0.05 — 0.1 is a sort 
of general rule. For the Salpeter case (soUd and dotted lines in 
Fig.|), Ar - 0.01 could be reached if assuming a larger star 
formation efficiency, /s « 0.1 — 0.2. 

However, before drawing any conclusion, it is necessary to 
perform another important check, related to the concomitant 
production of heavy metals, AZ. In fact, Pop-Ill metal enrich- 
ment should be constrained not to exceed observed upper lim- 
its, such as the lowest metallicities measured in Population-11 
Halo stars, or alternatively the metallicities measured in high- 
redshift, possibly little evolved, environments like the Lyman-a 
systems. 

From spectroscopic analyses of Pop-11 stars we can rea- 
sonably take Z « 10"^ as representative upper limit (see 
e.g. Norris et al. 1996). As for the high-redshift information 
the situation is more complex as reviewed by Pettini (1999, 
2000). At redshift z ^ 3, the estimated metallicity goes from ^ 
0.3-0.1Zq inLyman break galaxies (LBG), ~ O.I-O.OIZq in 
Damped Lyman-a systems (DLA), down to 0.01 — O.OOIZ© in 
Lyman-a forests. The metallicity values measured in LBGs and 
DLAs are relatively high and comparable with the metallicities 
measured in today's Pop-1 and Pop-11 stars, so that these sys- 
tems are likely to have already experienced significant chem- 
ical evolution. More plausibly, low-density Lyman-a forests 
may exhibit, instead, the genuine nucleosynthetic signature of 
Pop-Ill stars. 

Putting together all these pieces of information, in the 
present discussion we assume that the maximum allowed metal 



enrichment should be confined within 

10-3. 



10- 



Az„ 



From inspecting Fig. 10 



it turns out that for the Salpeter 
case (solid and dotted lines) the permitted interval of AZ„iax 
corresponds to a minor helium enrichment 10^^ ^ AF ^ 
few lO^"^, which is consistent with a quite low star formation 
efficiency, in the range lO^^*^ ^ /s < 10~^. This result agrees 
with previous findings by Ober et al. (1983), and Abia et al. 
(2001). 

On the other hand, in the case of the IMF peaked at M ~ 
1000 Mq (short-dashed line), AZ,„ax implies a significant he- 
lium enrichment, ISY > 0.01, and a much larger efficiency of 
star formation, /g > 0.1. An intermediate situation applies to 
the case of the extremely top-heavy IMF extending up to super- 
massive objects (long-dashed line). However, it is worth notic- 
ing that under the hypothesis that the first single burst of star 
formation produced just VMOs with M ^ 1000 Af© or super- 
massive objects, the primordial gas would have been enriched 
only with CNO elements (and possibly some a-elements). 
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without any contribution in heavier metals, like the iron-group 
elements. 

Finally, we remark that the helium-to-metal enrichment, 
AY/AZ, is independent of the star formation efficiency, but 
is affected by the individual stellar yields and the IMF features. 
According to our calculations, the two lower curves of Fig. |lO| 
correspond to a Ay/AZ ~ 0.15 — 0.6, respectively. These 
values are quite low compared to recent observational deter- 
minations based on HII regions and predictions of chemical 
evolution models, e.g. AF/AZ 1.9 ± 0.5 can be considered 
a representative value following Peimbert & Peimbert (2001). 
On the other hand, the highest curve is described by a huge 
helium-to-metal-enrichment ratio. Ay/ AZ w 900, because of 
the exceedingly-large helium stellar yield compared to that in 
the form of CNO elements. 

It is clear that this is an extreme situation and large room 
is allowed for intermediate results in between the case of the 
Salpeter IMF and that of a IMF peaked at very large masses. 
Ours is meant as an explorative test indicating that if the pri- 
mordial IMF was quite different from the classical one, the 
primordial chemical enrichment may have followed a different 
path from that one expected under standard prescriptions. 

4. Concluding remarks 

In this work we have discussed the evolutionary properties 
of zero-metallicity very massive objects (120 Mq < M < 
1000 Mq), on the basis of new calculations that extend the 
work presented in Paper I. Stellar isochrones (see Appendix 
A) for ages from 16 Gyr dow n to 10^ yr are available at the 
web-address 



http://pleiadi.pd astro.it 



In the attempt to estimate the possible effects produced 
by stellar winds from primordial VMOs, we adopt recent for- 
malisms that describe the role of radiation pressure and stellar 
rotation as mass-loss driving factors. The emerging picture is 
the following: 

- At extremely low metallicity (Z w 10^^), the mechanism 
of hne-radiation transfer should be scarcely efficient, ex- 
cept for very large masses, say ^ 75OAf0. We also find 
that, as long as the mass loss front penetrates into the chem- 
ical profile left by the convective core during the H-burning 
phase, the maximum CNO abundance exposed at the sur- 
face does not exceed sa 10~® — 10^*, which are typical 
values required by the CNO-cycle operating in stars with 
original metal-free composition. Such degree of chemical 
self-pollution is also too low to trigger efficient radiation- 
driven winds from zero-metallicity massive stars. Instead, 
much larger surface enrichment may be attained whenever 
the nuclear products of He-burning, like carbon and oxygen 
- potentially able to produce large mass-loss rates - were 
brought up to the surface by either the penetrating mass- 
loss front, or some dredge-up process. In our computations 
this occurs only for the 1000 AIq model. 

- Our calculations also indicate that rotating very massive 
stars may actually reach the critical condition defined as 
the ilF-limit, which should be likely accompanied by large 
mass-ejection rates (here taken as large as Merit = lO"'^ 



yr^ ^). However, the net impact on mass loss should be quite 
limited, given the extremely short time during which these 
critical regimes can be maintained. To this respect, we re- 
mind these results are obtained in the context of a simpli- 
fied description of stellar rotation, and the reader should 



consider the remarks expressed in Sect. 2.5 



- Finally we recall that, according to a recent analysis by 
Baraffe et al. (2001), another possible mass-loss driving 
mechanism, namely the pulsation instability - usually at 
work in VMO models with "normal" chemical composi- 
tion - has been found of modest potential efficiency at zero 
metallicity, at least for masses M ^ 500 Mq. 

By combining the supernova yields available in the litera- 
ture with our predicted wind contributions, we also evaluate, in 
a simple way, the chemical enrichment produced by a primeval 
population of (very) massive stars. It is found that a significant 
helium enrichment, AY ^ 0.01, may be reached by assum- 
ing that the primordial conditions in the Universe were such 
that only very massive stars could form, with typical masses of 
« 1000 A/0, or extending into the super-massive domain (up 
to w IO^Mq). 

In the former case the implied efficiency of gas-to-star con- 
version should be of the order of 10%, and the first metal en- 
richment should be in the form of CNO elements (no iron- 
group elements). On the other side, under the assumption of 
the standard Salpeter IMF, the resulting AY is indeed negligi- 
ble over the maximum allowed enrichment in metals, 10^^ ^ 
AZinax ^ 10^'^. The corresponding star formation efficiency 
should be also be quite small, « 0.01 — 0.1%. 

In any case, as already pointed out long ago by Bond et 
al. (1983) and Carr (1994), the interesting possibility arises 
that the first generation of stars could mask the true (some- 
what lower) primordial abundance of helium as predicted by 
the standard Big Bang Nucleosynthesis (e.g. Olive et al. 1997). 
The question and its implications have been recently addressed 
by Salvaterra & Ferrara (2002) on the base of the present re- 
sults. 



Appendix A: Isochirones 

Basing on the set of evolutionary tracks for VMOs presented in 
this work, isochrones for very young ages, from log t/yr = 6.6 
down to 4.0, have been constructed. In this way we comple- 
ment the set of isochrones for zero-metallicity stars already pre- 
sented in Paper I. In fact, for ages logt/yr < 6.6 the turn-off 
point corresponds to stars with M < 100 A/q, which defines 
the maximum initial stellar mass in the former set of tracks. 

An example of the extended isochrones in the H-R diagram 
is illustrated in Fig. A.l For ages 2 x 10^ ^ t/jr ^ 10^ the 
termination point marks the stage of central carbon ignition. 

Complete tables with the isochrones can be obtained 
upon request to the a uthors, or through the WWW site 



tittp://pleiadi. pd.astro.it They cover the complete age range 
from about 10"^ yr to 16 Gyr (4.0 < log(t/yr) < 10.25). 
Isochrones are provided at A log t = 0.05 intervals. We refer to 
Paper I (appendix B) and the quoted web-site for all the details 
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5 4.5 4 



log 

Fig.A.l. Theoretical isochrones in the HR diagram for the 

initial composition [Z = 0,Y = 0.23]. Age ranges from 
log(t/yr) = 4.0 to 6.6, at equally-spaced intervals of A log t = 
0.05. In aU isochrones, the main sequence extends down to 0.7 
Mq (log L/Lq ~ -0.5) 

concerning the description of the table format and layout, as 
well as the sources for magnitude and colour transformations. 
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